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Abstract 

The data on atmospheric neutrinos can be explained assuming the existence of 
oscihations between f^'s and a Hght sterile neutrino with mixing close to maximal, 
and 6m? ~ 3 x 10~^ eV^. This interpretation of the data is in potential conflict with 
the successes of big bang nucleosynthesis (BBN), since oscillations can result in a too 
large contribution of the sterile state to the energy density of the universe at the 
epoch of nucleosynthesis. The possibility to evade these cosmological constraints 
has been recently the object of some controversy. In this work we rediscuss this 
problem and find that the inclusion of a small mixing of the sterile state with Vr 
can result in the generation of a large lepton asymmetry that strongly suppress the 
i^/i ^ oscillations eliminating the possible conflict with BBN bounds. In this 
scheme the mass of the tau neutrino must be larger than few eV's and is compatible 
with cosmological bounds. Our calculations is performed using a Pauli-Boltzmann 
method. In this approach it is also possible to develop analytic calculations that 
allow physical insight in the processes considered and give support to the numerical 
results. 



1 Introduction 

The data on atmospheric neutrinos [|l], H have shown that muon neutrinos and antineutri- 
nos 'disappear' oscillating with a large mixing parameter sin^ 2^ ~ 1 and a squared mass 
difference |(5m^| ~ 3 x 10~^ eV^ into an 'invisible' state. Experimentally this state can 
be either a r or a sterile neutrino. The existence of light sterile neutrino states is also 
suggested by the results on solar neutrinos and of the LSND experiment ||^ that give 
independent hints for neutrino oscillations. 

If the observed oscillations are into a sterile state, this could have significant cosmo- 
logical implications, and in fact it has been argued ^, |^ that oscillations with these 
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parameters are already excluded by the present bounds obtained from primordial nu- 
cleosynthesis. Sterile neutrinos not mixed with ordinary neutrinos would have a small 
average density for temperature T < 100 MeV, and therefore have a negligible influence 
during the epoch of nucleosynthesis. However oscillations between sterile and standard 
neutrinos can be a source of sterile particles, and indeed if the oscillation have sufficiently 
large amplitude and are sufficiently rapid, the sterile states can be brought in thermal 
equilibrium with ordinary matter and give a significant contribution to the energy density. 
The success of Big Bang Nucleosynthesis (BBN) in predicting the relative abudances of 
primordial nuclear species can be used to put limits on the energy density during the 
epoch of nucleosynthesis or on the effective number of light particles present during this 
epoch. For example recently Olive, Steigman and Walker p| quote a (two sigma) limit 
< 3.3, assuming the so-called low-deuterium option. Such a limit can be translated 
into an allowed and a forbidden region for the parameters sin^ 29 and that describe the 
oscillations of a sterile neutrino with an ordinary one. Presently there is some controversy 
on the value of the limit for A''^''^; if the limit is relaxed to Nf^ < 4, as someone conserva- 
tively supports , obviously no region in the neutrino oscillation parameter space can be 
excluded, however it appears possible that a sterile neutrino in thermal equilibrium with 
ordinary matter could be excluded by BBN, once the controversial datum on primordial 
deuterium abundance from high redshift quasar absorbers will be definitely settled [1^ . 

Assume now that Nf^ ^ 4 is excluded by the considerations on nucleosynthesis: does 
this also exclude the interpretation of the atmospheric neutrino data in terms of z/^ 
oscillations ? We will argue here that this is not the case providing an explicit simple 
model where the cosmological bound is respected. 

In this work we will assume that a sterile neutrino state is mixed with the z/^ with 
oscillation parameters compatible with the fits to the atmospheric neutrino data, and 
we will compute the contribution to the energy density at the time of nucleosynthesis 
due to Us produced by oscillations. As already found by several authors ||, H, 0], if the 
sterile state is only mixed with the z/^, for the oscillation parameters suggested by the 
atmospheric z/ data, to a good approximation, it is indeed put in thermal equilibrium 
with the rest of the plasma, and the effective number of neutrinos during nucleosynthesis 
becomes Nf^ ~ 4, in possible conflict with the data on primordial nuclear abundances. It 
is however remarkable that the density of the sterile state depends critically on the details 
of its mixing with other ordinary neutrinos. In particular we will show that the inclusion 
of a very small mixing of the sterile state with the r neutrino (sin^ 26'^-^ ^ 10^^-10^^) can 
reduce dramatically the energy density of the sterile state to a neglibly small value if the 
neutrino that is prevalently z/,- has a mass m,y > 1 eV. This effect has been demonstrated 
for the flrst time by Foot and Volkas [|1^ . Qualitatively the mechanism that is operating 



is that Ur ^ oscillations develop earlier in time and generate a net r-lepton number 
in the plasma that acts to block the z/^ ^ Vs oscillations. 

This work is organized as follows. In section 2 we describe the effective potential 
for neutrinos in the hot and dense medium of the early universe and discuss how the 
oscillation parameters for mixing between standard and ordinary neutrinos are modifled 
by the presence of matter. In section 3 we discuss the method we have used to compute the 
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evolution with time of the populations of standard and sterile neutrinos. We have chosen 
in this work not to use the full machinery of the quantum kinetic equations but a simpler 
(Pauli-Boltzmann) approach, that is approximately valid when the standard neutrino 
mean free path is much shorter than the oscillation length. In section 4 we will discuss 
the evolution of the sterile neutrino density in the presence of a two-flavor mixing with a 
standard neutrino, considering the lepton number per unit comoving volume as constant. 
In section 5 we rediscuss the two-flavor mixing problem taking into account the fact 
that since neutrinos and antineutrinos in matter can oscillate with different amplitudes, 
the oscillations can generate a large lepton number in the medium, that acts back on 
the oscillations modifying the effective oscillation parameters. In section 6 we discuss 
the flavor evolution of a system of two standard (z/^ and z/^) and one sterile neutrino 
states, and determine the region of oscillation parameters that is in agreement with the 
^ij. ^ '^s interpretation of the atmospheric neutrino data, and also results in a small 
contribution of sterile neutrinos and anti-neutrinos to the energy density of the universe 
during nucleosynthesis. 



2 Effective neutrino mixing in the early universe 

We will study the cosmological consequences of assuming the existence of a two-state 
mixing between two different flavour neutrinos, z/q, and z/^, with a vacuum mixing angle 
6q (defined in the interval < ^ ^/4) and a squared mass difference = ml — m^. 
This means that the flavor eigenstates are a linear combination of the mass eigenstates 
^a, (with well defined masses nia and rrih) 

\ua) = cos 9o\i^a) + slu 9o\ub) , Wfs) = COS 9q \ z/fe) - slu 6'o | Z/q) • (1) 

In this way, for ultrarelativistic neutrinos of momentum p propagating in vacuum, the 
probability that a state |z^(t)) being a pure |z/q,) at an initial time is measured as a 
\up) at a time Un + r is: 

Poil^a Typ,tin + r) = \ {Uf3\u{Un + t))\^ = siu^ 2^o " slu^ (^^^^ (2) 

where we defined the vacuum oscillation length as iolp) = 2p/5m?. If the neutrinos prop- 
agate in matter, the effect of coherent forward scattering of neutrinos with the particles 
of the medium results in a flavor dependent effective potential that modifies the mixing 
angle and the oscillation length |T^: 

= ^'l"' (4) 

^/S'^ + {C-Va{p)+Vp{p)f 
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where we introduced the adimensional effective potentials Va,i3{p) = {2p / 5rn?')Va,p and we 
defined s = sin 26^0, c = cos 26*0- Consequently, for neutrinos propagating in a homoge- 
neous medium of constant density and composition one has that: 

Po{Va J^P, tin + r) > Pm(z/a ^ Z//3, tj„ + t) = siu^ 29m " siu^ (^2£~Jp)^ 

If the condition c = Va{p) — vp{p) is satisfied, then a resonance occurs and mixing is 



maximal . The resonant amplification of the oscillations was first studied in the solar 
environment were oscillations develop coherently and MSW conversion takes place. In 
the following we will study the mixing of a standard neutrino (z/q, with a = e, fi, r) and a 
strictly sterile neutrino, f ^ = f = in the previous expressions. The effective potential of 
a standard neutrino depends on its momentum and on the temperature and composition 
of the medium. 

We will study the effect of neutrino oscillations during the time when the universe has 
a temperature T < 150 MeV. This is the range of temperatures in which the oscillations 
can have effects on the products of BBN. The description of neutrino oscillations for higher 
temperatures and during the quark-hadron phase transition is a difficult problem, in any 
case the oscillations at high temperature are suppressed by matter effects as it will be 
discussed in the following. Note also that the population of sterile neutrinos produced by 
oscillations at high temperature would be afterwards diluted by photon production due 
to the decrease in the number of degrees of freedom. 

We will assume for simplicity that the comoving density of photons is constant, and 
that the density of muons in the medium is negligible. These assumptions are not exactly 
correct in the range of temperatures we are considering, but are good approximations. 
The ratio of the comoving density of muons and electrons N^{T)/Nf,{T) is 0.8 for T = 
~ 105 MeV, and drops to 0.1 for T ~ 26 MeV, as the entropy of the muons is 
transferred to the other particles of the medium, including the photons. Considering that 
the statistical weight of muons is about 20% when they are fully present, we can expect 
that the accuracy of our analysis is 0.02-0.2 for 25 < T/MeV < 150. On the other hand 
for T > 0.5 MeV ^ rrie the variation in the electron particle number can be neglected 
with a much better accuracy. 

The effective potential for a standard neutrino of flavor a can be written [|14|, |T5| as 
the sum of two terms 

t;„(p,T,L(")) = v^{p,T,L^'^^)+vl{p,T) (6) 
dm? 5m? 

The effective potential for antineutrinos is obtained changing the sign of the term, 
while the term remains unchanged. Note that in a medium with a net charge number 
L*^") neutrinos and anti-neutrinos have different oscillation properties. In equation we 
have defined y = p/T, and L^"^ is a linear combination of the charge asymmetries (per 
photon) of the medium that is relevant for the oscillations of standard neutrinos of flavor 
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a. For muon neutrinos we have 

L(^)=2L,^ + L,^+L,^-^5„, (8) 

where Qx = {ux — nx)/n^ with Q = L ot B. In equation (H) we have assumed that the 
electric charge of the medium vanishes (Le = Bp), and that all particles are in thermal 
equilibrium with the same temperature T. The asymmetry relevant for Ur ^ Vs oscilla- 
tions can be obtained from equation (|]) simply interchanging the labels jj, and r; for L^^'^ 
one has also to include the contribution of charged current interactions with electrons 
adding a term Le- The constants ai and are: 

aL = 2V2GF ^ 8.014 X 10-^2 MeV-2, (9) 

and 

be = b^\^l + 2j^j ^ 2.17 X 10-2° MeV-^ (11) 

Looking at equations (^ and (0) we can see that if the charge asymmetry of the 
medium is small, the mixing of neutrinos in matter is strongly suppressed: oc {T* /T)^"^ 
for temperatures much larger than a characteristic temperature 



T*(5m^)=(^) . (12) 



Numerically we have T* = 23.4 (18.9) \5m?'{eY'^)\6 MeV, for the oscillations of z/^'s and 
i/^'s (t'e's) with sterile neutrinos. 

2.1 The resonance condition 

In the following it will be very important to consider which neutrinos (or anti-neutrinos) 
have an effective mixing in matter that is maximal. This is verified for c — Va = 0, or 
more explicitely for 

5m? 5m? 

where the upper (lower) sign applies for neutrinos (anti-neutrinos). This equation is a 
relation between the temperature T, the asymmetry L^"") and the neutrino adimensional 
momentum y = p/T. As an illustration in fig. |l] we show for some representative values 
of the oscillation parameters the set of points in the plane {T,L^°'^) for which neutrinos 
and antineutrinos with adimensional momentum y = p/T = 2.2 (the maximum of the 
equilibrium Fermi distribution) have maximal mixing. 

Except for the special case of maximal mixing in vacuum, = 1, it is important to 
distinguish the cases of 5m'^ > and 5m'^ < 0. A positive (negative) 5m'^ means that 
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the mass eigenstate that is prevalently a sterile neutrino is the heavier (lighter) one. If 
< (and only in this case) both neutrinos and antineutrinos can have maximal 
mixing in an exactly symmetric medium (L^"^ = 0). For neutrinos with adimensional 
momentum y this happens at a temperature: 

TLiy) - (^) ' = T*i5m^) y-'s (14) 

Note that when the temperature decreases, neutrinos with higher and higher momentum 
y oc are resonant. More in general (always in the case < 0), given a certain 
value of the temperature T, for a non vanishing value of the asymmetry L^"^) , neutrinos 
and anti-neutrinos have maximal mixing for different adimensional momenta. One has: 



For a neutral medium the two values coincide, for L^"^ > one has y^^ > y^^, and 
the opposite for L'^"^ < 0. For large asymmetries of the medium the expressions for the 
resonant values of the adimensional momenta take the asymptotic forms: 



■'a 



and 

res _ c\Sm'\ (^j. 

For L^"^ positive yl^g^ applies to neutrinos and yl^^ to anti-neutrinos and vice versa for 
negative. 

In the case of Sm^ > 0, the situation of maximal mixing can be present only for 
neutrinos or antineutrinos (but not for both), and this is possible only for a sufficiently 
large charge asymmetry of the medium. The resonance condition is verified for neutrinos 
(antineutrinos) only for L^"^ > L* (L*-"'' < —L*) where 

L* = hl^KI^ (18) 

If this condition is satisfied the neutrinos (anti-neutrinos) have resonant mixing for two 
values of the momentum given by 



(the ± gives the two solutions for either neutrinos or anti-neutrinos). 

In the case of = 1, the sign of is not well defined, and the resonant condition 
can be obtained as the limit c — > from both the solutions found for positive or negative 
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dm?. Effective maximal mixing in matter is obviously present for both neutrinos and 
antineutrinos for y = 0, for any conditions of the medium, since in this situation all 
matter effects vanish; a second resonant value of the adimensional momentum y'^'^^ = 
{(iL/bfj,){\L^°'^\/T'^) exists for neutrinos (for L'-"^ > 0) or anti-neutrinos (for L'-"^ < 0). 

Some examples of the dependence on the momentum y of the mixing parameter in 
matter are shown in fig. |l^, to be discussed in section p. It is important to note that for 
a given condition of the medium (fixing the temperature and quantum numbers of the 
plasma) the width of the resonance, that is the momentum range for which neutrinos or 
antineutrinos have a large mixing is in general very narrow {6y/y ^ 1) and only a small 
fraction of all neutrino momenta can oscillate efficiently. Note also that the neutrinos 
have a non-negligible population only for values of y not too far from ?/peak — 2.2, and 
therefore if the resonant momentum is v''*''' n -C 0.1 or y'^'^^ n ^ 10 the resonance has 
little practical importance. 

Because of the narrowness of the resonant regions it is essential to study the flavor 
evolution of the neutrino ensembles considering in detail the momentum distribution of 
the particles. The approximation of considering monochromatic neutrino populations 
with a single momentum p = (p) does not give a sufficiently accurate description of the 
problem. 

3 Kinetic equations for the sterile and active neu- 
trino distributions 



A rigorous description of the flavor evolution of the ensemble of neutrinos in the early 
universe has to take into account at the same time the quantum coeherent effects of 
oscillations and the effects of interactions with the particles of the dense medium. This 
can be done rigorously making use of a density matrix formalism |T^. The interactions 
result in a loss of quantum coherence for the neutrinos [O, |18[ . It is possible to write the 



reduced density matrix for a system of two mixed neutrino flavors in the form 



(20) 



where aj (j = 1, 2, 3) are the Pauli matrices and P = (Pi, P2, P3) is a vector. The diagonal 
entries of p are the ratios of the classical distribution function (/^.^ (p) for z/q, and fu, (p) for 
Us) with the zero chemical potential Fermi-Dirac distribution feq{p) = 1/[1 + exp(p/T)], 
so that: 



Po{p) 



(21) 



Starting from flrst principles it is possible to derive [0 the following set of equations for 
P and Pq, dubbed quantum kinetic equations (QKE): 



di 



P(p,t) = yip,t)AP{p,t)-D{p,t)PT{p,t) + j^Poip,t)z 



(22) 
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^Pofet) = R{p,t) (23) 

where: 

Pt = P.x + PyV (24) 

V = + VyV + V;5 (25) 

V, = ^s, Vy = 0, V, = -^[c-Va{p,t)] (26) 

Analogous relations hold for antineutrinos. 

The use of local momenta p = p (R/R^) {R is the scale factor, R^ its reference value 
at a temperature T*) in the evolution equations (^2]) allows to take into account the 
expansion of the universe avoiding the explicit presence of a partial derivative term in the 
left hand sides (see for example 0]). As we are dealing with constant photon number, 
R/R^ = T^/T and thus we can use the convenient adimensional momentum y ^ p/T. 

The first term in the right hand side of equation (p2|) describes the coherent evolution 
of mixed neutrinos, while the second term describes the decoherence effect of interactions 



through a damping function D{y) that has the simple expression pl|: 



D{y) = , (27) 

where 

TM = GlKT'y (28) 

is the total collision rate for an active neutrino z/^ with momentum p = yT and ka ~ 0.92 
(1.27) for a = fj,,T {a = e). In equation (|28|) we have neglected Pauli blocking factors, 
and assumed that the charge asymmetry of the medium is negligible. Correction terms 



for the presence of asymmetries in the neutrino lepton number are given in |^ and are 
small. With this assumption the damping functions for neutrinos and anti-neutrinos are 
equal. It is also convenient to introduce an adimensional damping function: 

d^{y) ^ D{y) = ^^1^. (29) 



The third term {dPo/dt)z in equation (p^) takes into account the variation of the to- 
tal number of standard plus sterile neutrinos. This happens because, while z/q, oscillate 
producing z/^, the inelastic collisions refill the quantum z/^ states to preserve chemical equi- 
librium and this effect changes only the ^-component. The repopulation function R{y), 
driving the evolution of Pq, for T > 1 MeV has the form [ pT| : 



Riy) - rM < fo^^y^ > ■ (30) 

In equation (|30D , /f^ = 1/[1 + exp(?/ — C,a)] is the Fermi-Dirac distribution and = 
fia/T is the adimensional chemical potential; the approximation symbol means that the 
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expression is rigorously true only in the limit of complete thermal equilibrium for the non 
oscillating particles. When the average collision rate is much larger than the expansion 
rate of the universe H{T), the deviations from equilibrium are so small that the term 
{dPQ/dt)z can be neglected and equation (^2]) becomes a damped Block equation, well- 
known in studies of nuclear and electronic spin motion. An important thing to notice is 
that the initial condition of the system is necessarily P^- = 0,Pq = 1, because at high 
enough temperatures the damping parameter would be always so high to destroy any 
transversal component ( Turing or Zeno Paradox) jl^ and because the presence of sterile 
neutrinos is initially negligible. In fact a sterile neutrino decouples at very early times 
when T 3> 100 GeV and thus its density is highly diluted by the subsequent disappearance 
of degrees of freedom, especially at the quark-hadron phase transition. 

In this work, to study the flavor evolution of the ensemble of neutrinos, we will use a 
simpler set of equations that is given below: 



d 

di 



f,^{y^t) = T^s{y)-[hMt)-hMt)] 



dt 
d 



j^USy.t) = ^Uy) ■ [Us{y,t) - UAy.t)] + Y.^Ah 



(31) 

(32) 
(33) 

(34) 



where fll 



Ta{y) 



4 s2 + dl{y) + [c - vl{y) - v^{y)Y 



(35) 



^as{y) — 2 



1 + 



r f 



n -1 



As^ + dl{y) + [c-v]^{y) + vL 



(36) 



^0 and 



(for the explicit form in the right-hand side we have used the fact that = 
the definition ( pQ] ) for da)- The Cr[ ] in equations (p3D and (Q) are the collisional operators 
[pO|] for the different reactions involving the standard neutrinos and other particles in 
the medium. This set of equations can be rigorously derived from the quantum kinetic 



equations (QKE) under an appropriate set of approximations, as first discussed in |TT 



and shown in great detail in 



The equations have a simple probabilistic interpretation, and can be derived with 
heuristic considerations [^, [Tl| . Active and sterile neutrinos are described as two different 
particle species and oscillations work as a special inelastic process that can transform an 
active neutrino into a sterile neutrino and vice versa. The rate for the flavor transitions 
{a —>■ s OT s ^ a) can be derived considering oscillations that develop only between 
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interactions (or "measurement processes") that happen with an average time interval 
r = 2 /Vet (see Jill ^ more detailed discussion). From these considerations one obtains: 



r — r — '^'^^y^ I p ( + \\ — '^<^^y^ 2 

J- OS sa 2 \^rn\ya ^ ^sj f^in i T ) / j- 



p « \ 2 



(37) 



The solution of the flavor evolution equations is greatly simplifled assuming that the 
standard neutrinos remain in chemical equilibrium. The effect of the coUisional operators 
for sufficiently large temperature, i.e. for T > 3.5 (2.5) MeV for a = /x, r (a = e), is in 
fact to maintain, to a very good approximation, in thermal and chemical equilibrium the 
population of standard neutrinos. The elastic reactions redistribute the lepton number in 
a thermal way and the inelastic processes refill continously the quantum states depleted 
by the oscillations. 

The assumption of thermal equilibrium for the standard neutrinos does not determine 
completely their momentum distributions, because in general z/'s and I7's have different 
oscillation probabilities; this can generate an asymmetry in the lepton number L^^ (that 
cannot be cancelled by collisions) and a non vanishing chemical potential. The follow- 
ing relation connects the lepton number L^,^ to the adimensional chemical potentials for 
neutrinos and anti-neutrinos: 

L _^ f d'P ( ^ ^ \ 

which for ^a,^a ^ 1 becomes simply: 

Lu. = - + ^(Ct) (39) 

If the standard neutrinos are in thermal equilibrium one has = —^a and the active 
neutrino distributions are described only through L^^. Below the chemical decoupling 
temperature this is not true in general and one should take into account also the quantity 
^a + ^a'- this is important if one wants to describe the effects of electron chemical potentials 



on BBN p^. For the temperatures relevant for our discussion chemical equilibrium is 
verified, and therefore these effects are not important and we will make the approximation 
that a single quantity L,^^ fully determines the momentum distributions of both neutrinos 
an anti-neutrinos of flavor a. 

The simplifled framework we have discussed and that we will will use in this work gives 
a good description of the flavor evolution of neutrinos in the early universe only if some 
conditions are satisfled. A detailed discussion of these conditions and a formal derivations 
of the Pauli-Boltzmann equations from the full set of Quantum Kinetic Equations for the 



case of two-flavor oscillations of one active and one sterile neutrino can be found in ^1 
Here we will recall these conditions including some simple comments. There is a small 
angle condition, that is simply 

« 1 , (40) 
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or, when the effective mixing in matter is large (that is close to a resonance) the condition 



/ d 



a 



Each one of these conditions assures that all neutrinos in the early universe medium 
remain very close to a pure flavor (standard or sterile) eigenstate. This is obvious in 



the case of small effective mixing in matter; the condition ( ^Tf ) assures that frequent 
interactions interrupt the coherent flavor evolution, continuously reprojecting the neutrino 
states into flavor eigenstates. Therefore if the conditions (^O]) or ( ^l]) are satisfied the 
transverse components of the vectors P and P that describe coherent superposition of 
flavor eigenstates remain always small, the number of degrees of freedom is reduced, and 
the Pauli-Boltzmann equations become a good description of the evolution of neutrino 
populations. 

An additional requirement for the validity of the framework considered is that the 
relevant properties of the medium (such as the effective potential and the damping func- 
tion) do not change appreciably in the time £int between collisions, otherwise the diagonal 
terms of the density matrix have still an appreciable dependence on the derivatives of 
the off-diagonal terms. In ||2T| it is shown that the parameter variation is maximal at 
the resonance and that the (most restrictive) condition assumes the form D > QH. This 
condition implies that in order to be confident in the approximation the bulk of (anti- 
) neutrinos should enter the resonance at a temperature T^es ^ 3 MeV . If the charge 
number is negligible this implies that the Pauli-Boltzmann approach can be used only for 
|(5m^| > 10~^ eV^, see equation (0). For Tres ^ 3 MeV, full MSW conversion can occur 
at the resonance, an effect that the Pauli-Boltzmann approach neglects. In this situation 
major differences from a full quantum approach are expected. 

We note that even in presence of an appreciable transversal component of the vectors 
P and P approximately correct results can be obtained (provided that full MSW conver- 
sion does not occur) for the calculation of sterile neutrino contribution to the effective 
number of neutrinos plus antineutrinos. In fact the motion of the vectors around the V 
vector affects much more the difference (lepton asymmetry) than the sum. For the lepton 
asymmetry evolution the differences are more relevant and a classic approach misses the 
oscillatory behaviour at large angles [|T3], |23[ (of course below 3 MeV it also neglects 
MSW conversion effect on the lepton asymmetry evolution). Aiming at the evaluation 
of the effective neutrino plus antineutrino number it is meaningful to extend the classic 
approach also to high values of the mixing angle. To this purpose we will keep the term 
in the denominators of equations ( p5|) and (|36D . In proceeding directly from QKE, 



expressions for the rates that do not contain this term have been obtained. As we will see 
in next section, this would spoil the nice agreement of the results that we obtain using a 
Pauli-Boltzmann approach with those obtained using the QKE. 
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4 Sterile f production for constant L*^"^. 

In the previous section we have reduced the problem of the evolution of the ensemble of 
neutrinos in the early universe to the solution of a set of coupled differential equations. The 
flavor transition rates Vas and T^s depend on the lepton number L^^ that is determined by 
the oscillation themselves, and therefore we have a non-linear problem. In this section we 
will begin solving a simpler problem considering the asymmetry L^^ as constant. We will 
show that in the framework of our method, the solution of this problem can be expressed 
in a very simple form. In the next section we will discuss the complete problem and 
consider the lepton number dynamical variable. 

The first step for the calculation is to rearrange the equations in a more convenient 
form, replacing the distribution functions with the following quantities: 



J f 



eq 



fu. 

To 



(42) 



eq 



As we will deal with small chemical potentials (^^, .^q, <^ 1), we can perform the expansion: 



fuM=feq(^ + ^a 



ey + 1 



(43) 



Moreover, considering that the equilibrium distribution /^^ satisfies the Boltzmann equa- 
tion (see e.g. P 



Llfeq] ^ ^eqM - PH ^^fUv.t) ^ j/eqiy,t) 



it is easy to derive from equations 
Zs{y,t) (neglecting terms 0(^q)): 



IF) and 



0, (44) 
the following equations for Zs{y,t), 



d _ . , 



^as{y,t) 



^as{y,t) 



Uy.t) + iait)- 



(45) 



(46) 



ey + 1_ 

We are interested to calculate the contribution of sterile neutrinos to the energy density 
of the universe. For this purpose it is usual to define a sterile neutrino effective number: 

Pvs + Pp. 



2 flO 



(47) 



eq 



where p^^ = (7/8)(7r^/30)T^ is the energy density per degree of freedom of a fermion 
particle in thermal equilibrium with vanishing chemical potential. The results on BBN 
can be used to put constraints on this quantity at a temperaure T < 1 MeV. The sterile 
neutrino effective number can be expressed as: 



n: 



eflf 



120 

7^ 



dy 



r 

1 + ey 



4{y) 



(4J 
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where = {zs+Zs)/2. It is also interesting to consider the sterile neutrino charge number 



1 



2C(3) 



dyz 



r 



(49) 



where z~ = {z^ — Zs)/2. Thus we see that the interesting quantities to be calculated 
are zf and z~ , also because, when dealing with the coupled problem, these will be the 
quantities entering the equation for the a-neutrino lepton number. We can derive from 
equations (HS) and (HUf) the evolution equations for 2;+ and z~ , using = — 



dt 
dt 



r + r r 

^ as ^ as 4-1 , ^ as 

n [!-<] + — 



r — r 

^ as ^ as 



[1 



zT + 



r + r 

^ as ^ ^ c 



ey 



(50) 



(51) 



From equations (|35|) and ( p6D and using the definition = v{y,T) L^"'\ is easy to find 
the following expressions: 



r — r 

^ as J- c 



A 



v{c - vl) ■ 



r 4- r 

^ as ~ ■'- c 



2 A 



[s'+di+ic-vir + iv, 



.L\2-i 



(52) 
(53) 



with 



A^[s' + dl + {c-vl- v^J^] . [s^ + dl + {c-vl + v^^f] . (54) 

These equations allow to show that the terms explicitly containing C,a can be safely ne- 
glected: in the equation (|50|) the term containing is actually an C(^q) term, while in 
the equation (^TJ) the analogous term gives a small correction because [s'^ + d'^{y, T) + (c — 
v'^{y, T)Y + {v^Y] -C v{y, T) [c — v'^{y, T)] for the relevant values of y and in the range of 
temperatures we are considering. Thus we can neglect the terms in the equations for 
z, (|Sp and Zs (|§. 

It is convenient to adopt as independent variable in these equations the temperature 
instead of the time, using: 



^ = -HiT) T 
dt ^ ' 



Mp 



T 



where Mp is the Planck mass and 



an 



7r2 43 Stt ^ , , 

~ 5.44 . 

30 4 3 



One obtains in this way [|Tl 



^^s{y,T) 



^as{y,T) 

' H{T)T 



[l-Zs{y,T)] 



(55) 



(56) 



(57) 
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^Uy,T) = -^0^[i-Zs{y,T)] . (58) 

These equations are now simple enough to allow analytic solutions for the different pa- 
rameters involved (lepton number included). It is convenient to write the solution in the 
following form [p6| : 

z,{y,T) ^ = 1 - exp [-gsiy,T)] . (59) 

The function gs{y, T), assuming that the asymmetry is constant in time, can be calculated 
with a simple integral: 



oo 



9.(«,T) = J^ dr^-0n. (60) 

For Zg and one has simply to make the substitution T^s — > ^as- Using the definitions 
of Tas and H{T) one can write more explicitly: 

g,{y,T) = K^\6m^\-2 G^MT) , (61) 
with a similar expression for anti-neutrinos. In equation (pl| ) we defined 

r^2 K ]\/r 

= ^7 ^ ~ 898 (657) eV-^ (62) 
46q an 

for a = fi,T [a = e), and the function G^^{y, T), that depends on the parameters s^, Sm^ 
and : 

/•oo ,2 

G.Xy,T)= dt -2 (63) 



T/T{y,6nfi) ^2 + ^2 ^12 + 



c :p L(°)|5m2|-3 y-3 ± 



In equation (163|) the upper (lower) signs hold for 5m? > (5m^ < 0), and the constant 
do is given by 

(72 ^ 

(io = ^ 0.02 (0.008) for « = /i, r (a = e) . (64) 

The function G;7^ is obtained changing the sign of the term containing L*^°^ . In the lower 
limit of integration we have introduced the definition: 

ny,Sm')=(\p^y . (65) 



ba y'^ 

We are interested in the limit of small T, i.e. generally T < T{y, 5m?) (for typical y 
values). Note that for 5m^ < and in the absence of asymmetry, neutrinos and anti- 
neutrinos with adimensional momentum y have maximal mixing at T = cs T{y, 5m'^), see 
equation (|l4|). 



14 



4.1 Vanishing charge asymmetry of the medium 

It is interesting to consider first the case where the charge asymmetry of the medium 
is negligibly small. In this case the population of sterile neutrinos and anti-neutrinos 
created by the oscillations are identical. For L^°'^ = equation (|63|) simplifies to 

POO ^2 



where the positive (negative) sign refer to (5m^ > ((5m^ < 0), while neutrinos and anti- 
neutrinos have the same distribution. Note that all dependence on |5m^| and y is in 
the value of the lower limit of integration. For T — >■ the dependences on y and \Sm'^\ 
disappear and the function G depends only on the mixing parameter s^. 

poo ^2 

the ± sign again refers to the cases of positive or negative 5m^. 

In fig. 1^ we plot the functions -F-j-(s^). For maximal mixing the two curves have the 
same value {F±{1) ~ it/{6\/2) ~ 0.37), while for < 10~^ the two functions tend to 
values that differ by approximately two orders of magnitude (-F-(O) ~ 25.3, -F+(0) ~ 0.26). 
This is easily understood noting that for Sm? < and L^""^ = 0, both neutrinos and anti- 
neutrinos with adimensional momentum y go through a resonance for a brief time interval 
when the temperature is close to the value ce T{y, 5rn?)^ while in the case of positive 5rn? 
the resonance is absent. 

Most of the production of the sterile neutrinos for the case < happens during 
the short time when maximal mixing is produced. This is illustrated in fig. |^ that shows 
one example of the evolution of the sterile neutrino population for negative 5m?'. The 
fact that the final momentum distribution of the sterile neutrinos is proportional to a 
thermal distribution is the non trivial effect of an exact cancellation between the amount 
of time that each momentum spends at the resonance, and the degree of damping of the 
oscillations due to the interactions. 

Note that in the plane (s^, 5m?), a line of constant A^^''^ can be obtained from 



iV:f = l-exp 



For small s the corresponding line is: 



K^\5m'\^ s' F±{s') . (68) 



= log[l - <] . (69) 

In fig. I we show the curves in the plane (s^, |5m^|) that correspond to a constant value 
N^f that is proportional to the energy density of the sterile neutrinos at low temperature. 
It can be seen that the region of parameters that is allowed by the atmospheric neutrino 
data implies that the sterile neutrinos are fully thermalized (with Nf^ very close to 1) and 
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this is possibly incompatible with bounds from primordial nucleosynthesis. We note 
that this conclusion follows from the assumptions that the asymmetry of the medium is 
negligibly small and that the oscillations only happen between a sterile and a single active 
neutrino (and their antiparticles). 

Our analytic results improve those from the rough criterium Tas/H{T) < 1 ^, 



|2q , |g, 13] and can be compared directly with the numerical results from density matrix 
approach found in P, |^ (where momentum-dependence was not included, but it appears 
from our results that it does not play any role). In the case of interest for us, a = /x, r, 
the comparison gives a good agreement. It may be noted that for maximal mixing we 
have slightly too restrictive results since we neglected the effect of depletion of active 
neutrino numbers when transitions occur below the chemical decoupling temperature at 
about 3.5 Mev. For this same reason our results cannot be applied to the case a = e for 
5m? < 10^^. In this case the depletion effect must be taken into account due to the direct 
role of electron neutrinos in the neutron-proton interconversion rates and moreover the 
MSW transition plays a fundamental role. 

4.2 Non vanishing charge asymmetry 

We consider now the case where a non vanishing (constant) asymmetry is present in the 
medium. In this case the density of sterile neutrinos and anti-neutrinos produced by the 
oscillations could be generally different, and in turn also the net lepton number of standard 
neutrinos will not stay constant: it is reasonable to approximate it with a constant if it 
is varying sufficiently slowly, and in any case the complete, dynamical calculation will 
be considered in section ^ In fig. ^ we plot the effective number of sterile neutrinos at 
low temperature as a function of the (constant) asymmetry calculated assuming maximal 
mixing and for different values of \5rn? \ . It can be seen that for sufficient large values of L*^^^ 
the contribution of the sterile neutrinos to the energy density begins to drop following a 
power law A^^f oc [L(^')]-4 . From the fi gure we can see that with the oscillation parameters 
indicated by Super-Kamiokande the contribution of sterile neutrinos to the energy density 
would be sufficiently small if L'^'^^ > 10~^. 

It is not difficult to obtain analytically an understanding of the suppression of the 
oscillations due to the presence of the asymmetry L'^^\ For = 1, the denominator of 



the integrand in equation (63) can be rewritten: 

1 + + f[\L^^'^\/L±t^]^ (70) 

where the negative sign applies to the case of neutrinos for L^^^ > 0, to antineutrinos for 
L^^^ < 0, and we have defined 

2 

L = ^ \6m'\-^ y-s . (71) 

The results discussed in section [4.1| continue to hold if \L^^^\ <^ L; we consider here 
the opposite case, L^^^ 3> L (or —L^^^^> L). In this case the function Giy^(p^){y,T ^ 0) - 
see equation (^) - receives important contributions from two different regions: one is the 
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region of low temperatures, t ~ (L'^^-' / L) another is a high temperature region, where 
a resonance occurs for t = (L(^) / It happens that the first region (low temperature) 

gives the dominant contribution (and therefore the productions of sterile neutrinos and 
antineutrinos are similar), and this explains the behavior of the sterile neutrino number 
with increasing lepton asymmetry, visible in fig. ^ namely 

G'^:{y,T^ 0)1^,^^ ^ 0.43 {L^^^ L)-'/' . (72) 

The resonant high temperature contribution would give instead a behavior 

G^fiy,T^O)U^^±-[^^)~" , (73) 

that is therefore more suppressed at large L^^\ 

Note that for ~ 1 the effect of a charge asymmetry of the medium results in a 
suppression of the sterile neutrino and anti-neutrino population. The situation can be 
different for small. In this case the presence of a charge asymmetry can even result 
in an enhancement of the oscillations. This is illustrated in figs. ^ and |^. In fig. ^ we 
show a case of enhancement of the sterile production, present for Srn^ > and ^ 1 
(a similar behavior was also noticed for much smaller values of in ||29[). For these 
oscillation parameters, in a neutral medium there is never a situation of maximal mixing, 
however for a sufficiently large asymmetry neutrinos or antineutrinos with adimensional 
momentum y can pass through a resonance for two well defined values of the temperature 
(see fig. |T]). At the lowest T value the production of sterile neutrinos (or anti-neutrinos) 
can be significant, and in fact, for positive asymmetry, the production of neutrinos for 
positive becomes equal to the production of antineutrinos for negative 5m^, since in 
both cases we have resonances with almost identical characteristics. In fig. ^ the behavior 
for large L^^^ L shows a decrease with the same power law as in equation ([72|). 

For still smaller values of s^, it is possible to have an intermediate region of asym- 
metries, still larger than L, where the sterile neutrino production increases with the 
asymmetry. An example of this behavior is given in fig. |^. In fact, integrating with a 
saddle point approximation for L^^^L ^ 1, one finds for the resonant component (in the 
approximation c ~ 1): 

i_ 

2 

(74) 

Note that for very large L^^^ the oscillations are suppressed as before, however for a certain 
range of values of L^^^ determined by the value of one can have an enhancement. 

It should be clear from the arguments given above that the main contributions to the 
integrals come from rather restricted regions of temperature (particularly if a resonance is 
present) and therefore the requirement of a constant L^'^^ for the validity of these results 
can be slightly relaxed: it is in fact sufficient that the asymmetry is slowly varying in the 
relevant region of the integration. 
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5 Mixing of two neutrino flavors 



In this section we will study the flavor evolution of the neutrino populations in the case of 
two-flavor mixing of a standard and a sterile neutrino, considering the lepton asymmetry 
of the medium as a dynamical variable. We will conflrm the flnding of |2^, that for 



< one has in general the generation of a large lepton asymmetry Ly^ ; at the same 
time, summing over z/'s and z7's, the energy density of the sterile neutrinos generated by 
the oscillations remains similar to what has been obtained in the previous section using 
the (incorrect) assumption that the charge asymmetry of the medium remains small. 
The set of differential equations that we need to solve (|3T|) , (|3^) , (0), (|34D has been 



introduced in section 2. In section 3 we have recalled that the coUisional operators main- 
tain the active neutrinos in thermal equilibrium, therefore their momentum distributions 
are determined by the chemical potential iaiT) = —^a(T), or equivalently (see equation 
(|38D) by the lepton number L^^iT). Thus the flavor evolution is described by a set of 
only two equations, such as (^) and (|^) or (^Up and (^). We have however to take into 
account that the lepton number of the medium is evolving with time. The only source 
of variation of the lepton number is the presence of oscillations, therefore one has that 
Li,^ + remains constant and: 

L,^ (T) + const = -L,, (T) = ^ j dy (y, T) + (y , T)] . (75) 

In summary the problem of the flavor evolution of the neutrino population has been 
reduced to the solution of the set of two equations ( ^5]) and (|§), with the additional 
condition that the asymmetry L^j^ that enters in the expression of the transition rates Tas 
and determines the distributions of the standard neutrino and anti-neutrino population 



is obtained with the integral (|75|). 

It is instructive to consider an equation directly for L^^ that can be obtained combining 
the deflnition ([75| ) with equations (^5]) and (|4^). The result is: 



dL,^ 1 



dt 4C(3) 



j dy y" jl^{y) [(r,. - f„.)(4 - 4) + (r,, + f„,)(;^; - z-)] . (76) 



This equation allows to understand the evolution of the lepton number making use of some 
simple considerations. The second term tends to bring the neutrino population toward 
a situation in which the initial lepton asymmetry is equally shared between active and 
sterile neutrinos. The flrst term, for positive 5m^, has always sign opposite to L("), and 
also acts as a damping term, so that the asymmetry approaches a constant value —fi/2, 
with ff = L(") -2Ly^. 

On the other hand, if Sm^ is negative, the flrst term can have either the opposite or 
the same sign as L^"^ and may also be the source of a growth for the asymmetry. To have 
an appreciable growth we need that the transition rates for neutrinos and anti-neutrinos 
are different and that the total number of active neutrinos and sterile neutrinos are also 
different [Q. Therefore considering the oscillation of two active neutrino species that 
have approximately equal number densities, one expects only an effect of equipartition of 
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the initial asymmetry. A small effect would arise after chemical decoupling because of a 
small number density difference between the active neutrinos. This is created by electron- 
positron annihilations that reheat only the electron neutrino component |31]. Things can 
be much more interesting in the case of active-sterile neutrino mixing. 

In fig. ^ we show as a function of the temperature the energy density of the sterile 
neutrinos, always calculated assuming the existence of z/^(,-) oscillations with max- 

imal mixing, positive 5m'^ = 3 x 10~^ eV^ and different initial values of the asymmetry 
L*^") (0, 10""^, 10~^, 10~^ and 10~^). The thin curves show the results of a calculation 
performed considering L*^"^ constant in time: increasing L^°'^ the sterile neutrino produc- 
tion is more and more suppressed. The thick lines have been calculated assuming the 
initial values for the asymmetry, but taking into account the evolution with time of the 
lepton number. To a good approximation for any initial value L^""^ < 10""^ the initial 
lepton number is efficiently destroyed and the evolution of the effective number of sterile 
neutrinos is independent on the initial lepton number and similar to the L^""^ = case 
(three curves are coincident with the thick solid line). For L*^"-* > 10~^ the lepton number 
remains approximately constant and the results are very similar to the case in which the 
evolution of the lepton number with temperature is neglected. 

In fig. H we show the evolution with temperature of the lepton number L,^^ (a = /i, r) 
calculated with a numerical integration of the flavor evolution equations, assuming the 
existence of Ug ^/^(r) oscillations with a negative (5m^, and some representative choices 
of the mixing angle. The evolution of the neutrino populations is started at T = 150 MeV 
with Lf°^ = 10"^° and r/ = 5 X 10"^^ 

One can distinguish several phases in the evolution of L^^: (i) initially the asymmetry 
j^ia) jg erased and L,^^ is brought to the value —fi/2, (ii) then at a critical temperature 
Tc the asymmetry starts an exponential growth, (iii) after a short time the exponential 
growth is transformed into a power law growth, and (iv) finally the asymmetry is 'frozen' 
and remains at a constant value. 

Some insight on the onset of the instability can be obtained with the following qualita- 
tive discussion, where we will assume for simplicity that the oscillating flavors are and 
Us- In this case, as we have discussed in section 2, the resonant momenta for neutrinos and 
anti-neutrinos are ordered as yl^^ > y^^ for L*^"^) > or yl^^^ < y^^ for L*^"^) < (the values 
are equal in a neutral medium). Let us now consider the evolution of the asymmetry 
L^'^^ with temperature starting from an arbitrary initial value. The values of the resonant 
momenta, initially very low, increase with time yl'^^ ~ y-^^ ~ ^/c[T*{6m'^)/T)]^. When 
the resonant values become of order unity the oscillations begin to have effect and start 
to produce sterile neutrinos and anti-neutrinos. For a certain time, when the position 
of the resonant momenta is below the peak of the Fermi-Dirac distribution of the stan- 
dard neutrinos, no lepton asymmetry is produced and in fact any small initial non-zero 
value of L*^"^^ is erased. This can be easily understood observing that if in this situation 
(when both resonant values of the adimensional momentum are below the peak of the 
distribution) a small positive (negative) L^'^^ develops, the resonance for the neutrinos 
will move to higher (lower) values of y, with respect to the position of the resonance for 
anti-neutrinos; therefore the neutrinos will oscillate more (less) efficiently, and since the 
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sterile particles are much less abundant that the standard ones, this will result in the net 
disappearance of z/^'s (z7r's) and in the generation of a negative (positive) Ly^ that will 
tend to cancel the initial fluctuation of the charge number L^^\ In summary an initial 
value of the charge asymmetry will be erased and L^'^^ will take a value close to zero. 
This mechanism of destruction of an initial charge asymmetry cannot work if the initial 
asymmetry is too large and becomes inefficient for |L(")| > 10"^ pi. 



The argument outlined above to justify the stability of the condition L^^^ = obviously 
loses its validity when the resonant momenta are above the peak of the Fermi-Dirac 
distribution, at |/ ~ 2.2. More precisely when the temperature drops below a critical 
value T = Tc such that 

yl^:^yL^^^^c(^) ^2.2 (77) 



the situation becomes unstable. A small positive fluctuation of L^'^^ results in a larger 
resonant value for the neutrinos so that the oscillations of the anti-neutrinos become 
more efficient, more V^s disappear and the asymmetry L^,^ increases even more. The 
same argument shows that an initial negative fluctuation in L^^^ after Tc, would grow to 
a much larger (in absolute value) negative asymmetry. 

At flrst sight it appears that we have a genuine case of instability and that the sign of 
the lepton number will be determined only by random fluctuations. This is in fact not the 
case, at least in the theoretical framework we are considering. It is possible to demonstrate 
that asymptotically (for T — > 0) the sign of the lepton asymmetry will be the same as the 
sign of fj = Ly^ + Li, —Bn/2, which is the charge asymmetry of the medium excluding the 



contributions of the i/^'s and that in the present discussion remains constant |rT|. Note 
that it can be expected that fj is positive. This can be derived from the assumptions that 
(i) the electric charge of the universe vanishes (Lg = Bp), (ii) the net quantum number 
B — Loi the universe also vanishes, (iii) protons and neutrons have equal number densities 
(for T > 3 MeV) and therefore Bp = Bn, and (iv) the total lepton number of the standard 
neutrinos before the development of the oscillations is equally divided among all flavors 
{Ly^ = Ly^ = Ly^). From these assumptions one can derive f/ = -Bp/6, that is a positive 
number of order 10^^*^. At the critical temperature, where the term proportional to Ly^ in 



the right hand side of equation (|76D vanishes, and the sterile neutrino asymmetry is still 
negligible, the sign of the time derivative of Ly^ is the same as the sign of fj\ therefore, 
once the term driving the instability starts to dominate, \Ly^\ will continue to grow, and 
the sign will be the same as the sign of fj. We will come back to this point and briefly 
reconsider the limitations of the present approach in the concuding section. 

In flg. and |Tl| we show the evolution of the distributions of sterile neutrinos and 
antineutrinos with temperature. We can again distinguish different regimes. In a flrst 
phase, for T larger than a critical temperature T^, the initial small asymmetry of the 
medium is erased, the resonant momentum for neutrinos and anti-neutrinos grows {y'^'^^ oc 
T^^) producing sterile neutrinos with higher and higher momentum. During this phase (as 
discussed in section 4) the sterile neutrino population, for y < y^^^(T), is proportional to 
a thermal distribution, that is Zs{y) and Zs{y) are constants up to a maximum y < y^'^^(T) 
(compare with flg. H). 
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When the temperature drops below the critical temperature and the production of a 
positive Ly^ starts, the resonant momentum for anti-neutrino oscillations is pushed back 
to lower values producing additional sterile anti-neutrinos. In fig. 10 it appears clearly 
that the critical temperature corresponds to the situation where the resonant momentum 
was ?/^^;(-^) ~ 2.2. 

Fig. O is calculated for the same value of 5m?, but for a larger mixing parameter 
s^. In this case the critical temperature is lower and corresponds to a situation where 
the resonant momenta for neutrinos and anti-neutrinos is larger, ■0'''^% ^ ~ 7.1. An un- 
derstanding of this "delay" of the onset of the instability for large can be obtained 
inspecting equation (ffGl). It can be seen that the term that controls the growth of the 
asymmetry is {Vas — fas) (-2^ — z^). This term is reduced if the sterile neutrinos develop 
a sizeable population, as it is the case for a sufficiently large s^. We can see in fig. |Tl| 
that the sterile neutrinos produced in the early phase (T > Tc) are in fact a significant 
fraction of the thermal equilibrium population. 

The estimate of the energy density in sterile neutrinos and anti-neutrinos produced by 
oscillations for T — 0, is modified with respect to the estimates developed in the previous 
section, because of the effects of the evolution of the lepton asymmetry that determines 
the neutrino effective potential. It is important to note that the difference between the 
estimates (with and without the dynamical evolution of Ly^) is small for those oscillation 
parameters that predict an appreciable energy density, A*"^^ ^0.3. This is illustrated in 
fig. |T2| where the two thin solid lines in the plane (s^, 5ir?) correspond to N^^ = 0.1 
and 0.6 in a complete calculation of the fiavor evolution, while the two thin dot-dashed 
lines are calculated assuming a constant negligible L^,^. Note that there is some differ- 
ence in the results for A^^^ = 0.1 (so that in the complete calculation the region having 
A^^^ < 0.1 is larger), but the two lines corresponding to A^^^ = 0.6 are essentially undis- 
tinguishable. This can be understood qualitatively (see also the discussion in [|ll|, |33[|): 
when the production of sterile neutrino is sufficiently strong, the lepton number genera- 
tion is delayed (the critical temperature becomes lower), and the subsequent production 
of sterile neutrinos in the regime with a not negligible lepton number occurs far from the 



peak of the Fermi distribution, and therefore gives a small contribution, see fig. Ill 



Thus we can conclude that in the simplest framework of two-fiavor active-sterile neu- 
trino mixing a detailed treatment of the lepton number generation does not essentially 
modify the region of oscillation parameters that the constraints of BBN can exclude. As 
we will discuss in the next section, this is not true any more if we consider a more general 
framework with more elaborate forms of mixing between neutrino fiavors. 



6 Three neutrinos mixing scenario 



In this section we will consider the scenario proposed by Foot and Volkas ||rT| where a 
sterile neutrino is mixed with both i/^ and Ur- As in 1 11] we will not consider the most 



general case of three neutrino mixing but we will limit our considerations to a situation 



where we can neglect the 



oscillations and the descriptions of all transitions can 
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be well approximated by a pair of two-flavor oscillations: ^ Ug and Ur ^ Ug. In 
the scenario we are considering two neutrino eigenstates are nearly degenerate with a 
squared mass difference of order Smj^ ~ 3 x 10~^ eV^, and are both superpositions with 
approximately equal weights of z/^ and Us and only a small component of Ur] in this way 
the z/^ oscillations can describe the atmospheric neutrino data. The third mass 

eigenstate is nearly a pure Ur- 

As we have seen before, if the Uj. does not participate in the oscillations, the ^ Ug 
transitions result in an energy density in sterile neutrinos not compatible with suggested 
bounds from big bang nucleosynthesis, unless an unnaturally high lepton number L^^^ is 
present at the start and persists down to temperatures T < 10 MeV. The inclusion of a 
small mixing of the r neutrinos deeply modifies this conclusion. The oscillations between 
Uj. and Ug can produce a significant L^^ and this in turn can suppress the production of 
sterile neutrinos in the z/^ ^ Vg oscillations. 

For this mechanism to work it is necessary to require |5m^| > |5m^| (so that the 
Pg oscillations start earlier), 5m1 < (the mass eigenstate dominantly coupled to 
Vr is the heaviest) and s^. sufficiently small (to assure that the lepton number generation 
can actually occur). Moreover the constraint (|69D must be verified, otherwise a sterile 
neutrino overproduction would occur already during the early Vg oscillations. These 
conditions are not sufficient by themselves to prevent a subsequent destruction of the 
lepton number by the z/^j Vg oscillations and as a consequence N^^ ~ 1. In the following 
the conditions on the oscillation parameters to avoid this conclusion are discussed. 

We indicate with Sq,, Cq, 5rn?^ (a = /x, r) the mixing parameters for the z/q, Vg 
oscillations. The set of equations that describe the flavor evolution for the proposed 
scenario can be easily obtained with a simple generalization of the case of two flavor 
mixing discussed in section ^, describing at the same time the dynamical evolution of 
the two lepton numbers Ly^ and Ly^. As before the active (anti-) neutrinos are assumed 
to be in thermal equilibrium, namely, neglecting terms of second order in the chemical 
potential, z+ = z+ ~ 1, z'^f^^^ ~ 12 C(3) Lj,^j^j/[7r^ (1 + e^^)]. Thus we can write: 



dt 4C(3) 
dLy„ 1 



j dy y' f^^iy) { [T^g - T^g] (l - 4) + [r.. + T.,] {z; - z;)} (78) 

dy y' fliy) { [r,s - T,g] (i - zj) + [T,g + T,g] [z- - z-) } (79) 



dt 4C(3) 

j^zg{y, t) = [r„ + r^,] [1 - zgiy, t)] (80) 

j^Uy,t) = [r„ + r^s] [i -^.(y,t)] (8i) 

Note that the evolutions of the two asymmetries Ly^ and Ly^ are coupled for two reasons, 
the first being simply that Ug oscillate both in and in z/^, the second is that, also when the 
population of sterile neutrinos is small and the transitions Ug u^^r) can be neglected, 
the asymmetry of one type of flavor (for example L^^) is relevant for the evolution of 
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the other asymmetry (in this case Ly ) because it modifies the effective potential in the 



medium of the neutrino {v^) 



We have numerically integrated this set of equations calculating as a function of tem- 
perature the distributions of sterile neutrinos and the asymmetries of z/^'s and z/^'s. The 



results show some very interesting and not entirely intuitive features. In fig. |T2| we show 
curves of constant energy density in sterile neutrinos for T — >■ as a function of the 
oscillation parameters Sm^ and assuming = 1 and a fixed value of Sm^^. The most 
remarkable feature is that for a given value of 6m?^ the plane (s^, 5m^) is divided in 
two regions. In one (the region below the thick lines) the sterile neutrinos become fully 
thermalized, while in the other (above these lines) the energy density in sterile neutrinos 
drops to a much lower value. This value is almost undistinguishable from what is obtained 
in a calculation with only Vr ^ oscillations (the ^ Ug oscillations being completely 
damped). In fact, the lines of constant A^^''^ in fig. |12] drawn as thin lines are exactly the 
same that were obtained and discussed in section ^ using a two-neutrino scenario. 

The boundary in parameter space between the region where the sterile neutrino is 
thermalized (resulting in N^^ ~ 1) and the region where the effect of oscillations 
is essentially negligible is given to a first approximation by the relation 

4 (Smlf ^ 28 {6mlf (82) 

(see the dashed curves in fig. |T2|). From the figure, we read that the minimum value 
allowed for the tau-neutrino mass in order to comply the BBN bounds is ~ 1.4 eV (3 eV) 
if SttiI = 10~^ (3.2 X 10~^) eV^, with a mixing angle ~ 10~^. 



In fig. we show the evolution with temperature of the sterile neutrino effective num- 
ber considering values for the neutrino oscillation parameters that are very close to each 
other but belong to the two different regions where the sterile neutrinos are or are not 
thermalized. The thick lines correspond to a point immediately below the corresponding 
boundary in fig. O the sterile neutrino production is momentarily interrupted when the 



Ur asymmetry is generated at a temperature of about 30 MeV, but, since the initial rapid 
increase of L^^^ is stopped and then reversed, it starts again and reaches the thermal equi- 
librium at a temperature ~ 10 MeV. The thin lines, corresponding to a point immediately 
above the boundary in fig. 0, describe a situation in which L^^^ continues to increase, 
and the z/^ oscillations are damped and unable to generate further sterile neutrinos. 

In the following we will present some qualitative arguments and analytic estimates to 
illustrate these numerical results. 



6.1 Qualitative discussion 

In fig. |l^ we show an illustration of the neutrino mixing parameters in matter at different 
times (or temperatures). This particular example was calculated for the set of oscillation 
parameters: = 10^^ eV^, = 1, (5m^ = —10 eV^, = 10^^. In each panel the four 
different curves describe the effective mixing parameter for the four transitions between 
standard and sterile (anti-)neutrinos: s'j^ .^, sf^ ,^, s^,^, and s^-^; the different panels refer 
to different temperatures, from top to bottom the temperatures and asymmetries in the 
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five panels correspond to the situations (a) : T = 58.1 MeV, L^^ = 0, (b) : T = 38.7 MeV, 
L^^ = 0, (c) : T = 34.9 MeV, L^^ = 8.9 x 10^^ (d) : T = 32.2 MeV, L^^ = 1.7 x 10-^ (e) : 
T = 29.1 MeV, Li^^ = 3.2 x 10^^ and in all cases Lj,^ ^ 0. The behaviour of the evolution 
of the mixing parameters is easy to understand qualitatively. At high temperatures the 
medium is neutral and the oscillation parameters in matter for z/'s and I^'s are equal. For 
vanishing charge asymmetry of the medium, the mixing parameter in matter for z/^ ^ 
oscillations is 



2 —2 
s = « 
/i,m Ml™ 



T 



y~' (83) 



and the mixing is strongly suppressed by matter effects. For the Vj- ^ oscillations, the 
important feature is the fact that at a momentum 

C = C = (^) ^ T-' = V^r (^)' (84) 

the mixing for neutrinos and antineutrinos is maximal. With decreasing temperature 
the suppression of the z/^ ^ oscillations becomes weaker while the position of the 
resonance for the z/^ <-> Ug oscillations grows as T~^. When yl'^^ approaches the peak of 
the Fermi distribution, the instability discussed previously sets in and an asymmetry L^y^ 
becomes to develop. The positions of the resonances for the tau-sterile oscillations of 



neutrinos and anti-neutrinos are given by equation ([TsD for a=T. If a positive L^^ starts 
to be generated the resonance for the neutrinos moves to higher momenta y, where the 
neutrino population is smaller, therefore Vr oscillate more effectively, and the asymmetry 
grows faster. 

An important consequence of the generated L,^^ is to modify the oscillation parameters 
of the z/^ ^ Us oscillations. The presence of an L^^^ asymmetry depresses the oscillations 
for most values of the neutrino momentum, but also induces maximal mixing for neutrinos 
or antineutrinos (depending on the sign of L*^'^^) in a narrow range of momenta. For 
positive L(/^) ~ the resonance is present for neutrinos at a momentum 

C-^^ (85) 

and it gives rise to the generation of a negative L^,^, that also acts on the oscillation 
parameters for both type of oscillations. 

It appears clearly now that two different scenarios can develop. In the first case the 
generation of L,y^ is never sufficiently rapid to cancel the growth of L,y^ , and the combined 
asymmetry L(^) ~ 2L. + L,^ keeps growing. The resonant peak for the muon neutrino 
oscillations rapidly (yl^^ oc L'^^^T'^) moves to very high values of y where it becomes 
irrelevant because no neutrinos have such a high momentum, and the z/^ ^ z/^ oscillations 
remain suppressed by the generated L^^^ (this is the situation envisaged in fig. TZ). In 



the second possible outcome, the oscillations <-> Ug proceed so rapidly that at a certain 
time the growth of L^^"^ is stopped and then reversed. Qualitatively it is understandable 
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that in the space of oscillation parameters there will be a sharp boundary between regions 
where the two different scenarios develop. For example let us keep fixed Smf^, 8m^ and 
= 1, and consider the evolution of the neutrino populations for different s^. In a certain 
interval of the oscillations of tau-neutrinos proceed sufficiently rapidly to generate a 
net L'^^^ and suppress the muon-neutrino oscillations. Decreasing the rates for the 



transitions, equations (p5| , |3^) , become weaker and weaker and the generation of 
slower. At a critical value the Vr ^ Vg oscillations do not proceed sufficiently rapidly 
to cancel the effect of the <-> Ug oscillations and the growth of L^^^ can be stopped and 
reversed. 

Qualitatively it should be clear that for a set of parameters along the boundary between 
the two regions the conditions (temperature and asymmetries of the medium) when the 
derivative dL'^^^ / dt vanishes and the growth of the asymmetry is reversed are such that the 
resonance for the Vg transitions is close to the maximum of the Fermi distribution 

(jj^es ^ i/pg^j^ = 2.2), where it affects the largest possible neutrino population. 



6.2 Analytic estimate 

To obtain an estimate of the boundary that separates the regions in the space of the 
neutrino oscillation parameters where the sterile neutrino production is suppressed, we 
will make the assumption that the reversal of the asymmetry L'^^^ happens when both the 
resonances for Ug and Vj. <-> Vg oscillations, that are the source of a negative dL^^ 

and a positive dL^^ are most effective, that is in conditions where both resonances are at 
a value y ~ y^cak — 2.2 (this corresponds to the situation of panel (d) in fig. |T^. We will 
therefore impose the condition: 



dt dt 



= (86) 



This algorithm is motivated by our numerical results that indicate that the reversal hap- 
pens for values of the temperature and the asymmetries L^,^ and L^y^ such that this 
condition is approximately satisfied. The algorithm can also be understood qualitatively 
observing that if there is no reversal of the asymmetry L^''^ the resonance momentum yl^^ 
takes all values from zero to infinity and therefore it will take the value ypeak] the effect of 
the z/^ Ug oscillations (relative to the Ur ^ 'i^s oscillations) is on the other hand most 
important when y]^^ ~ y^eak- During the early phase of the generation of L^^ one also has 
y^^ < ?/peak, at least for sufficiently small (as illustrated in an example in fig. p!0| ). 

The expressions for the derivatives of the lepton asymmetries can be approximated, 
keeping only the dominant terms, in the following form: 



dL^^ 



dt 4C(3) 
dL„.. 1 



dt 4C(3) 



dyy' f(y) [F., - F^] (87) 



dy f(y) [T,g - T,g] (88) 
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Neglecting furthermore the oscillations of I7^'s that are non resonant and the oscilla- 
tions of the z/t-'s that are happening at large values of and calculating the integrals with 
the saddle point approximation since they are dominated by the resonance contributions, 
equation (^61) can be rewritten as: 



+ 4 



12 



-) 



„ {Gpkr 

4 + - 



2 u \2rpl2„A 



(„. _ bf^T 



-, -1 



1/peak 



|W| 



Vpeak \ 



-, -1 



LM J |5m2| 



0(89) 



For (GlKfT^^yl^JlSmll^ > and (G 



2 k ^2 



T 



12 



?/peak/l'^ 



2 |2 



^ 1 the dependence 



on Vpeak drops away (showing that the choice of that particular value was not critical) 
and the condition for the boundary takes the simple form 



2|2 



2 2 + 



\6m 



2 |2 



(90) 



In this way we have reproduced the functional form (s^ \Sm^ 
boundary between the regions in parameter space where the i/^ 



^'^ ~ const |(5m^p) of the 
^ Us oscillations develop 



or not (see fig. 0). The numerical results can be fitted with the constant approximately 
equal to 28, and with the analytic estimate sketched above we obtain a reasonable approx- 
imation for it. In fact, neglecting the small term L^^ — Bp/2, one has: L(^) = 2U + U 
and L^/^) = 2U + U ; since Li, is negative we have L^'^^ L^'^^ > 2. 



It is also possible to understand why the thick lines in fig. ^ are not straight for 
relatively large s^, and bend upwards. This is due to the fact that in this situation the 
resonant momentum y^^ at the critical temperature is larger than ypeak, and therefore the 
approximation of equal y"^^^ and y^'^^ fails. The constraint becomes more restrictive because 
the tau-sterile antineutrino oscillations are less effective, and a larger 6m'^ is needed to 
compensate this effect. 



6.3 Validity of the approximations 

It is appropriate at this point to comment again on the validity of the approximation 
scheme used in this work. As discussed in section |^, the Pauli-Boltzmann system of 
equations describes correctly the flavor evolution of the neutrino ensemble only if the 
effective mixing in matter is negligibly small or (for large mixing) if iint/^m "C 1, that is 
if the interaction length is much shorter that the oscillation length and the oscillations 
cannot develop coherently. We can now check a posteriori that these conditions are 
indeed well satisfied for neutrino oscillation parameters in the 'allowed region' indicated 
in fig. 12, during the crucial epoch when the development of the i/^ oscillations is 

decided; therefore the use of the approximation is valid, and the determination of the 
allowed region is correct. 
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As discussed before the crucial moment in the flavor evolution of the neutrino popula- 
tions is around a temperature T* ~ Cr\Sm'^\'s /br determined by the 6ml inass difference. 
Around this temperature — > and V^- Vg transitions begin to proceed with slightly 
different rates, and a positive L^,^ asymmetry starts to be generated. If the transitions 
i^fi do not proceed sufficiently rapidly to overcompensate the growth of L^^ when it 
is still very small, a large L^^"* ~ Ly^ develops that supresses the effective v^-Vs mixing 
for all lower temperatures. On the contrary, if the — > Vs oscillations are sufficiently 
effective to reverse the growth of L^^\ the medium does not develop a large asymmetry, 
and at lower temperatures (T < \5m?^\a jh^) the oscillations can proceed without being 
suppressed by matter effects, resulting in a large population of sterile neutrinos at the 
nucleosynthesis epoch. 

We can now investigate the validity of our framework for T ~ T^. For this temperature, 
the effective mixing in matter is very small (also in the case of the ^ Ug oscillations 
where = 1), except for narrow regions around the resonant momenta. Therefore to 
check the validity of our framework we have to check that iint/^m is indeed small for the 
neutrinos at the resonances. Using the definitions of ii^t and im we can write explicitely 
this condition as: 



\6ml\ 



< 1 (91) 



For the z/^ ^ Vs oscillations, for small asymmetry of the medium, (as it is the case 
for T ~ T*) the resonant condition is: hrT^ y'^ /\5ml\ ^ Cr, and substituting in (|9lD one 
obtains the condition 

^ « ^ 0.04 (92) 

Cr Ot 

In the allowed region we have s,- ^ 3 x 10"'^, and this condition is always satisfied. 

For the ^ Ug oscillations the mixing in matter is negligibly small for all momenta y 
except for small regions around y = 0, and in case of a non symmetric medium around a 
second (non vanishing) value, (B^). The first region is not significant because for T ~ T^, 



it is narrow (the mixing is large only for values y < y\Sm'j^/6m'^\) and for the values of the 

mass differences we are considering it contains a neglibly small population of neutrinos. 
For the second region, substituting T ~ the condition (0) can be rewritten as: 



1 

7^2 



6m 



drni 



< 0.04 (93) 



The boundary of the allowed region is approximately described by the relation (^), 
substituting |5m^/5m^| in ( P^D one obtains 

^ < 0.21 (94) 

y 

This condition is satisfied (for all y with a significant neutrino population) along most 
of the boundary line. The Pauli-Boltzmann approximation begins to be poor for Sr 
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large, that is only at the 'lower' corner of the allowed region where one can expect small 
deviations from the results obtained solving the detailed Quantum Kinetic Equations. 

We note that if the oscillation parameters are outside the allowed regions, and the 
medium of the expanding universe remains quasi-symmetric, at sufficiently low temper- 
ature the matter effects become negligible and the z/^ Ug oscillations can proceed 
coherently with maximal mixing. In this situation the Pauli-Boltzmann approximation 
is not formally valid, however in this case the final result is equal populations of sterile 
and standard neutrinos, independently of the method used. 

The validity of the Pauli-Boltzmann approach is also very well confirmed by the 
calculation of the boundary of the allowed region performed using the Quantum Kinetic 
Equations by Foot |2^, whose results are in good agreement with what we obtain here. 
The only small deviations (the bound is slightly more stringent than ours) are present in 
the lowest part, where we expect our approximation to be less accurate. 

We stress again that the advantage of simpler equations, allowing the use of analytical 
arguments to check the numerical results, is hard to overestimate in a complicated problem 
such as the one we are considering here. 



7 Discussion of the results and conclusions 



Our results are in good agreement with the numerical results obtained in [|Tl|, |25[ , except 
in the region of the highest acceptable s^. The curves of [|ll|, calculated in a Pauli- 



Boltzmann approach, do not show an upper bending as ours do, while the results obtained 
in 



25|] using QKE are slightly more restrictive than our curves. 
Our results do not agree with the semi-analytical calculations presented in |^ , that 
obtain a much more stringent lower limit on the tau-neutrino mass {rriu^ > 15 eV), that 
would give too much energy density in the form of hot dark matter and problems for early 
structure formations. The authors of [0, ^ claim that very small values of yl'^^ are most 
important to determine the fate of the lepton number and that an insufficient resolution 
in the momenta may well miss this contribution and thus yield wrong results. We have 



accurately tested our numerical results; for example, the curves in figure |I3| describing 
the behaviour of the asymmetry L^^^ remain the same, if we change the minimal value of 
y used to calculate the sterile neutrino distribution down to 10~^. This strengthens our 
result, that the values of yl^^ that can better destroy the lepton number are around i/peak 
and not the small ones. We performed the same test in several points in parameter space. 



that essentially cover the thick curves in figure 12. Moreover we note that for small values 



of yl^^ (more precisely, for yl^^ < y^^^ y\6m'j^\/\6m'^\ ~10 '^ — 10 ^ according to the region 

of parameters considered), the sterile neutrino production from z/^ — Ug transitions has 
a dominant non-resonant behaviour (and cannot therefore be evaluated in the resonant 
approximation), which is rather smooth and does not require a resolution particularly 
high to be estimated correctly. 

Furthermore, we do not agree with the condition imposed in |33, R^f to determine the 



range of parameters that allow the survival of the lepton number generated. In those 
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papers the amount of lepton number destroyed by z/^ ^ Ug oscillations for any yl°^ during 
the time needed to cross the resonance is compared with the total value of the lepton 
number L^^^ produced up to that moment. In our opinion, the comparison should be 
done instead with the amount of lepton number created by z/^ ^ Ug oscillations during 
the same time interval (i.e. making a comparison of the two rates as we did). Actually 
this amount is comparable with the total lepton number created before and in fact it is 
even more, since the lepton number is exponentially growing. 

The same authors have also recently found another argument that should restrict the 



allowed region |3g] . The QKE show a chaotic behaviour for the sign of the lepton number 
that is generated at the critical temperature ||22|, ^ and this would lead to the formation 
of lepton domains. Neutrino crossing the borders of these domains would pass through a 
further resonance, with an additive contribution to sterile neutrino production. We think 
that even though on qualitative grounds this observation sets forth an important issue, 
from a quantitative point of view only a full momentum dependent approach can give a 
definitive answer about the chaotic behaviour (moreover, the present results of momentum 
independent calculations ^ are not in good agreement with each other). In ||2^ such 
an approach is used and the chaotic behaviour has also been found, but only for large 
enough mixing angles. A complete analysis of the region in parameter space where chaos 
occurs is still missing however, and we think that the effect of lepton domain formation 
on the allowed region may be more accurately determined. 

In conclusion, a momentum dependent Pauli-Boltzmann approach for active-sterile 
neutrino oscillations in the Early Universe allows to derive, directly from the equations, 
useful analytical results for the sterile neutrino production that support the numerical 
ones. An analytical procedure provides a deeper physical insight and allows a more gen- 
eral picture of the possible solutions that sometimes pure numerical calculations hide (the 
generation of a relevant lepton number missed in early numerical results is a good exam- 
ple). In this way we have been able to derive analytically in section ^jthe constraints from 
BBN on mixing parameters with the assumption of negligible lepton numbers, showing 
that they essentially reproduce the results that can be obtained from QKE for z/^(^) ^ Ug 
oscillations. When lepton numbers are present their effect is for most cases a suppression 
of sterile neutrino production, but we showed that it exists a regime where an enhance- 
ment is possible, even if in a parameter region not interesting for the present BBN bound. 
We were also able to give a full spectral description of the production, both in the case 
of negligible lepton numbers and when lepton numbers must be considered. This allowed 
to specify better why in the case of a two neutrino mixing the regions allowed by BBN 
do not change significantly with or without a generation of lepton number. The most 
important result has been obtained about the possibility that the generation of lepton 
number through tau-sterile neutrino oscillations can suppress sterile neutrino production 
from mu-sterile oscillations for the values of mixing parameters that could describe the 
Super-Kamiokande results. We confirmed numerically previous results [|ll|, that show 
that this is possible with a lower limit on the tau neutrino mass (of a few eV) compatible 
with present cosmological observations. We could provide an analytical support to these 
results and show that one should not worry to obtain an exceedingly high resolution in 
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momentum space. This settles in our opinion the recent controversy on the minimum Ut 
mass allowed. 

We did not consider the more general situation where lepton domains can be created. 
We think that clear results about this issue are necessary, but still missing. Anyway it is 
surely important to have established results in the simpler homogeneous case. 

Acknowledgments 

We would like to thank E.Kh. Akhmedov, R. Foot and R.R. Volkas for useful discussions. 



30 



References 

[1] Super-Kamiokande collaboration, Y. Fukuda et ai, Phys. Rev. Lett. 81, 1562 (1998); 
K. Scliolberg (for the SK collaboration), |hep-ex/9905016 . 



[2] Kamiokande collaboration, Y. Fukuda et ai, Phys. Lett. B 335, 237 (1994); 
1MB collaboration, R. Becker-Szendy et a/., Phys. Rev. D 46, 3720 (1992); 
Soudan 2 Collaboration, W.W.M. Alhson et ai, Phys. Lett. B 449, 137 (1999); 
MACRO collaboration, M. Ambrosio et ai, Phys. Lett. B 434, 451 (1998); 
F. Ronga (for the MACRO collaboration), |hep-ex/ 98 10008 . 

[3] See for example J.N. Bahcall, P.I. Krastev and A.Yu. Smirnov, Phys. Rev. D 58, 
096016 (1998) and references therein. 

[4] LSND Collaboration, C. Athanassopoulos et al, Phys.Rev. C 54, 2685 (1996); 
Phys. Rev. Lett. 81, 1774 (1998); Phys. Rev. C 58, 2489, (1998). 

[5] R. Barbieri and A.D. Dolgov, Phys.Lett. B 237, 464 (1990); Nucl.Phys. B 349, 743 
(1991). 

[6] K. Enqvist, K. Kainulainen and M.J. Thomson, Nucl.Phys. B 373, 498 (1992). 

[7] X. Shi, D.N. Schramm and B.D. Fields, Phys.Rev. D 48, 2568 (1993). 

[8] K.A. Olive, G. Steigman and T.P. Walker, |astro-ph/ 9905320] , to be published in 
Physics Reports (1999). 

[9] E. Lisi, S. Sarkar and F.L. Villante, Phys.Rev. D 59, 123520 (1999). 



[10] D. Kirkman, D. Tytler, S. Buries, D. Lubin and J.M. O'Meara, |astro-ph/9907T28 
and references therein. 

[11] R. Foot and R.R. Volkas, Phys.Rev. D 55, 5147 (1997). 

[12] L. Wolfenstein, Phys.Rev. D 17, 2369 (1978). 

[13] S.P. Mikheyev and A.Yu. Smirnov, Nuovo Cimento C 9, 1 (1986). 

[14] D. Notzold and G. Raffelt, Nucl. Phys. B 307, 924 (1988). 

[15] K. Enqvist, K. Kainulainen and J. Maalampi, Nucl.Phys. B 349, 754 (1991). 

[16] A.D. Dolgov, Sov.J.Nucl.Phys. 33, 700 (1981). 

[17] R.A. Harris and L. Stodolsky, Phys.Lett. B 116, 464 (1982). 

[18] L. Stodolsky, Phys.Rev. D 36, 2273 (1987). 

[19] B.H.J. McKellar and M.J. Thomson, Phys.Rev. D 49, 2710 (1994). 



31 



[20] J. Bernstein, "Kinetic Theory in the Expanding Universe", Cambridge University 
Press, Cambridge, 1988. 



[21 
[22 
[23 
[24 
[25 
[26 
[27 
[28 
[29 
[30 

[31 
[32 
[33 
[34; 
[35 
[36 
[37; 



N.F. Bell, R.R. Volkas and Y.Y.Y. Wong, Phys.Rev. D 59, 113001 (1999). 
X. Shi, Phys.Rev. D 54, 2753 (1996). 

R. Foot, M.J. Thomson and R.R. Volkas, Phys.Rev. D 53, 5349 (1996). 

R. Foot and R.R. Volkas, Phys.Rev. D 56, 6653 (1997). 

R. Foot, Astropart.Phys. 10, 253 (1999). 

J.M. Cline, Phys.Rev.Lett. 68, 3137 (1992). 

D. Fargion and M.G. Shepkin, Phys.Lett. B 146, 46 (1984). 

K. Kainulainen, Phys.Lett. B 244, 191 (1990). 

D.P. Kirilova and M.V. Chizhov, Nucl.Phys. B 534, 447 (1998). 

M.Yu. Khlopov and S.T. Petcov, Phys.Lett. B 99, 117 (1981); Erratum: Phys.Lett. 
B 100, 520 (1981). 

P. Langacker, S.T. Petcov, G. Steigman and S. Toshev, Nucl.Phys. B 282, 589 (1987). 
R. Foot and R.R. Volkas, Phys.Rev.Lett. 75, 4350 (1995). 
R. Foot and R.R. Volkas, Astrop.Phys. 7, 283 (1997). 
X. Shi and CM. Fuller, Phys.Rev. D 59, 063006 (1999). 
X. Shi and CM. Fuller, |astro-ph/ 98 12232. 



X. Shi and CM. Fuller, |astro-ph/990404T 



K. Enqvist, K. Kainulainen and A. Sorri, |hep-ph/9906452 . 



32 




5 10 15 20 25 30 



T (MeV) 



Figure 1: Curves in the plane (T,L*^"^) that correspond to maximal mixing for neutrino 
and anti-neutrinos with y — p/T — 2.2. All lines are calculated for \5'w?\ — 1 eV^ and 
« 1. 
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Figure 2: Plot of the functions (dashed hue) and F_{s'^) (sohd hue) as a function 

of the mixing parameter. 
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Figure 3: Evolution of the sterile neutrino distribution assuming the presence of ^ 
^fiir) oscillations with = — 1 cV^ and = 10^^. The flavor evolution is followed 
from an initial temperature Tj = 150 MeV, where we assumed that sterile neutrinos are 
absent. We plot the ratio of momentum distribution of i^^'s and Vg^s (equal in this case) 
with a thermal equilibrium distribution with zero chemical potential. The thick solid line 
shows the distribution at T = 0.1 McV, the thin dot-dashed, long-dashed and dashed 
lines show the distributions for higher values of T. At each temperature the production 
of Us occurs only in a small interval of momenta where oscillations are resonant. As an 
illustration the dotted curve shows the region where sterile neutrinos are produced for T 
just below 150 MeV. The direction of the arrows indicates how the resonant momentum 
is changing with time. 
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Figure 4: Lines in the plane (s^,|5m^|) that correspond to a fixed value of the sterile 
neutrino effective number A^^^: 0.1, 0.3, 0.5, 0.7, 0.9 (from bottom up). The curves are 
calculated assuming that the asymmetry of the medium remains negligibly small. The 
solid (dashed) lines correspond to < (5m^ > 0). The thick solid line shows the 
allowed region indicated by the Super-Kamiokande experimental results. 
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Figure 5: Plot of the effective number of sterile neutrinos A^,^^(eff) as a function of L^^^ 
calculated assuming that the sterile neutrino is maximally mixed with a muon (or a tau) 
neutrino. The curves are calculated for |5m^| = 10~^, 10^^, 10^"' and 10^^ eV^ (from 
right to left). In the calculation the asymmetry L^^^ is fixed and it is not considered as a 
dynamical variable. 
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Figure 6: Plot of the effective number of sterile neutrinos A^j,^(eff) as a function of L'^^^ 
calculated assuming that the sterile neutrino is mixed with a muon (or a tau) neutrino with 
— 10~^ and \6m'^\ — 1 eV^. The thick (thin) lines correspond to Sm^ > {5m^ < 0). 
The solid (dashed) curves correspond to neutrinos (antineutrinos) for L^'*) > and vice 
versa for L^'*) < 0. In the calculation the asymmetry L^'*) is fixed and it is not considered 
as a dynamical variable. 
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Figure 7: Evolution with temperature of the density N^f and asymmetry L,^^ of sterile 
neutrinos, calculated assuming the existence of Ug ^ ^nir) oscillations with Sm^ = — 1 eV^ 
and = 10^^, and a time independent value for the charge asymmetry of the medium. 
The thick (thin) lines describe the evolution of iV^^ (-^i^J- The dotted lines correspond to 
L^"^ = 0, the dashed to L^^*) = 10"^ and the dot-dashed ones to = 10"^. For L^'^) = 0, 
i/'s and 17's oscillate in the same way and no asymmetry is generated. For non negligible 
values of L^^^ the sterile antineutrino production can be enhanced, and an asymmetry is 
generated. 
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Figure 8: Evolution of the sterile neutrino energy density calculated assuming the exis- 
tence of Vg oscillations with maximal mixing and Sw^ = 3 x 10~^ eV^. The solid, 

dotted, dashed, long dashed and dot-dashed lines correspond to an initial value = 0, 
10"^ 10-^ 10"^ and 10"^. The thin lines are calculated considering L^"-* as a a constant, 
the thick lines taking into account the dynamical evolution of L'^"\ 
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Figure 9: Evolution with temperature of the lepton number L^^ {a = /i, r) (thick hnes) 
calculated assuming the existence of <-> oscillations with the indicated oscillation 
parameters. The evolution of the neutrino populations is started at T = 150 McV with 
= 10^^*^ and fj = 5 x 10^^^. The thin lines are the resonance curves for the indicated 
values of momenta (dashed for neutrinos and dot-dashed for antineutrinos) . 
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Figure 10: Evolution of the sterile neutrino distribution with temperature for z/^(^) 
oscillations when lepton number is generated. The oscillation parameters are = 
—1000 eV^ and = 10"^ and correspond to the thick long-dashed line in figure ^. The 
different curves correspond to increasing temperature. Above the critical temperature 
(Tc ~ 60 MeV) a standard production with zero lepton number occurs, up to y^'^^ ~ 2.2. At 
lower temperature, an enhancement regime with constant y^'^^ ~ 0.35 is clearly apparent. 
The production occurring after lepton number generation gives a dominant contribution 
in this case, but the total production is anyway too small to be relevant for BBN bounds. 
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y=p/T 

Figure 11: Evolution of the sterile neutrino distribution with temperature for ^ti(T) 
oscillations when lepton number is generated. The oscillation parameters are 5m? = 
—1000 eV^ and = 10~^ and correspond to the thick dotted line in figure |^. The 
different curves correspond to increasing temperature. Above the critical temperature 
(Tc ~ 42 MeV) a standard production with zero lepton number occurs. It can be noticed 
that in this case y^'^^ ~ 7.1, since z+ < 1. The enhancement regime now is too short to be 
observed and the dominant contribution to the total production comes from the standard 
regime with zero lepton number. 
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5rn^). The thin sohd 

hnes correspond to A*"^^ = 0.1 and 0.6; the thin dot-dashed hnes correspond to the same 
values of A^^''^ and have been calculated in a two neutrino mixing framework with negligible 
L'^'^\ The three thick solid lines, calculated for the three values = 10~^, 10~^'^ and 
10~^ e"V^, correspond to the points where the energy density of sterile neutrinos passes 
quasi-discontinuosly from the value A*""^ = 1 (below the curve) to a much lower value 
(above the curve). The three dashed lines, corresponding to equation (^2]), give a first 
approximation description of the numerical results. 



Figure 12: Lines corresponding to constant A'^^^ in the plane (s:;. 
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Figure 13: Evolution the asymmetry L^^^^ (solid lines) and the sterile neutrino energy 
density N^^ (dashed lines) as a function of the temperature. The thin curves are calculated 
for the set oscillation parameters = 10"^'^ eV^, and = 10~^, and = 17.80 eV^, 
The thick curves are calculated for the same values of and but with a slightly 
smaller value = 17.40 eV. 
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Figure 14: Mixing parameters of i^'s and I^'s in matter for different conditions of the 
medium. The oscillation parameters are: = 10~^ eV^, = 1, = —10 eV^, 

= 10^^. The 4 panels refer to the situations: (a) T = 58.1 MeV, L^^ = 0, (b) 
T = 38.7 MeV, L^^ = 0, (c) T = 34.9 MeV, L^^ = 8.9 x 10"^ (d) T = 32.2 MeV, 
L^^ = 1.7 X 10-^ (e) T = 29.1 MeV, = 3.2 x 10"^; in all cases L^^ ~ 0. In each panel 
the thick (thin) curves refer to z/^ i^s {I'r ^ J^s) oscillations, the solid (dashed) lines to 
neutrinos (antineutrinos). The thin, smooth curve is a plot of the function feqiu)- 
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